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Abstract 

A non-equilibrium Green's function method is applied to model high-field quantum transport and 
electron-phonon resonances in semiconductor superlattices. The field-dependent density of states 
for elastic (impurity) scattering is found non-perturbatively in an approach which can be applied 
to both high and low electric fields. I — V curves, and specifically electron-phonon resonances, are 
calculated by treating the inelastic (LO phonon) scattering perturbatively. Calculations show how 
strong impurity scattering suppresses the electron-phonon resonance peaks in I — V curves, and 
their detailed sensitivity to the size, strength and concentration of impurities. 
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I. INTRODUCTION 



There has been tremendous interest in the electron transport properties of semiconductor 
superlattices (SLs) since the pioneering work of Esaki and Tsu 1,2 . Bragg scattering of carriers 
from the Brillouin zone boundaries of a crystal produces Bloch oscillations, or, equivalently, 
localized (Wannier-Stark) states within each unit cell period if the system is sufficiently clean. 
This regime is characterized by negative differential conductance (NDC) at sufficiently high 
fields, and occurs if the Bloch frequency Q = eFL/h (F is the electric field and L the unit 
cell period) is larger than the effective scattering rate l/r e //. The large unit cell of the 
superlattice can allow electrons in a SL to perform these Bloch oscillations before being 
scattere d 3 i 4 i 5 i 6 . 

If the electric potential drop across n unit cells corresponds to an optical phonon en- 
ergy, Bryksin and Firsov^ predicted another high-field effect: so-called electron-phonon 
resonances, where an LO phonon mediates a transition between localized Wannier-Stark 
states whose energy separation equals the phonon energy. These resonances manifest as 
peaks in the current-voltage relationship and produce a nonmonotonic current-voltage de- 
pendence. Experimental I — V characteristics of cubic ZnS films have been identified with 
such resonances^. Growing interest in the transport properties of semiconductor superlat- 
tices, e.g. GaAs/AlAs and InAs/GaSb, has resulted in several recent theoretical studies 
of these electron-phonon resonances in SL a 9 ^ 11 ' 12 ! 13 ? 14 . Unfortunately, to our knowledge, 
electron-phonon resonances have not been observed experimentally in SL transport within 
the NDC regime. 

Several choices are available to calculate SL transport in the high and low field regimes. 
Wacker and Jauho^ 5 - established the relationship in some cases between different transport 
calculation methods, including miniband conduction^, Wannier-Stark hoppingi 6 -, sequential 
tunnelling^, and the non-equilibrium Green's function (NGF) formalism^ and their respec- 
tive ranges of validity, and reached the conclusion that the NGF contains the other three as 
limiting cases. For example, the semiclassical Boltzmann transport equation (BTE), which 
can be solved to any desired degree of numerical accuracy by Monte Carlo simulations^ 9 -, 
cannot accurately describe transport under high-field (eFL ~ Tiuoq) condition al 20 . Use of 
the generalized Kadanoff-Baym(GKB) ansat z 18 ! 21 ' 22 is limited in validity to slowly varying 
fields and time dependence. 
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Rott, Linder and Dohler— employed the Wannier-Stark hopping transport technique 
to study the temperature dependence of electron-phonon resonances. Their method is only 
applicable to the high-field regime. Emin and Hart^ obtained a 5— function electron-phonon 
resonance in the limit of a very wide energy band A —>■ oo. Bryksin and Kleinert focused 
on narrow band(A — > 0) transport in a one- dimensional NGF model with 10 and without 
the GKB ansatz. The approximations made in Ref. lj| included relying on the periodic 
nature of the Green's functions in a large homogeneous field and discarding the Fourier 
components of the Dyson equations beyond lowest order. The impurities were modeled with 
5-function potentials (and hence only depended on the product of the impurity density and 
the scattering potentia l 11 ! 15 ), and the SL bandwidths considered were much smaller than 
optical phonon energies, in order to obtain analytic expressions for the field-dependent DOS 
and the transport. 

In our treatment here we use the Kadanoff-Baym-Keldysh NGF techniqu e-""'-^ to calculate 
the non-linear quantum transport. We solve for the field-dependent Green's functions equa- 
tions numerically, non-perturbatively including the elastic impurity scattering by directly 
inverting the NGF Dyson equations. These Green's functions, and the resulting density of 
states (DOS), depend on the impurity size, strength, and density separately. We then cal- 
culate the effect of inelastic (phonon) scattering as a perturbation. Specifically we leave the 
retarded Green's function obtained for elastic scattering unchanged, but alter the less-than 
Green's function. A related model, where impurity scattering was treated exactly within 
the resonant level model, and phonon scattering was treated perturbatively, has been intro- 
a,eed befe,* .» Ref. Q *e anabs, focused on comparison of the model results with 
those obtained from the GKB ansatz. In our treatment we consider more general impurity 
types, and use the solvability of the model to explore the influence of impurity nature on 
the electron-phonon resonances. The resulting I — V curves illuminate the nature of the 
electron-phonon resonances in a regime where elastic scattering dominates. 

Because we treat more of the problem numerically we are able to consider all significant 
Fourier components of the Dyson equation, general SL bandwidths and impurity potential 
size and strength. In this way we can explore the properties of these electron-phonon reso- 
nances for a wider range of parameters with confidence. Our approach therefore possesses 
several advantages over earlier ones: (i) it employs a more realistic form for impurity poten- 
tials that leads to better approximations for the DOS; (ii) other common approximations 
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are relaxed, including restrictions on the growth axis miniband width; (iii) it is applicable to 
both high and low field transport; and (iv) it does not use the generalized Kadanoff-Baym 
ansata 1 ^ to solve the Dyson equation. 

The outline of the paper is as follows. In Sec. |H1 we show the detailed theoretical and 
numerical methods used to calculate the retarded and lesser Green's functions, from which 
are extected the DOS and carrier distribution functions. In Sec. IIII[ the field dependence of 
the calculated drift velocity is discussed in detail. We focus on the electron-phonon resonance 
peaks and the influence of impurity scattering and establish a relationship between earlier 
results and ours. Our results are discussed and summarized in Sec. IIVI 



II. THEORY AND NUMERICAL METHOD 



A. High-field Green's function 



In order to study high-field effects, we must use an approach that is non-perturbative in 
the electric field. We consider the scalar potential gauge 0(r, t) = — r ■ F, where F is the 
time-independent uniform electric field. In the absence of scattering, the Dyson equation of 
the retarded Green's function in the k representation can be written as^ 3 - 



d 

ih— - e(k) + ieF • V k 
at 



G r Jk, t, k', t') = (27r) 3 <5(k - k')5(t - f). 



(2.1) 



We make use of a 14-band superlattice K ■ p calculation 24 to obtain the superlattice energy 
dispersion relations s{k). Only the lowest conduction band with growth-axis width A is 
considered in the present transport calculations. We assume the electric field is along the 
superlattice growth direction (noted as ||, and the in-plane directions are noted as _L). The 
K • p results are then fit to a dispersion relation of the form 
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where m* is the electron effective mass along the in-plane directions. Solving Eq. (j2.1j) . we 
get the retarded Green's function with no scattering 
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After taking the Fourier transformation (t — t' — > u 
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The corresponding gauge invariant spectral function is 



A(k,w) = z[G r (k,cu) -G a (k,cj)] 
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The spectral function satisfies the sum rule (H/2ti) J dujA(k, uj) = 1. The density of states is 
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The energy of the Wannier-Stark level n superlattice cells away is neFL. A level with index 
n can contribute to the DOS if the miniband width A is big enough, and its contribution is 
proportional to J 2 (— ^wz)- This result coincides with Wacker's results based on Wannier- 
Stark bands^. In the NGF approach of Bryksinii the terms with n > vanish because the 
approximation A — > is used. 



B. Elastic impurity scattering and G r 

For a single impurity at R, the Dyson equation for the retarded Green's function is 

G^(k,k» = <^(k,k» + fj^J ^^( k ^i^) 

xe - 4(qi - q2 ).R y(qi; qa> w )G^(qa, k', to). (2.7) 

For elastic impurity scattering, we have V r (qi,q2,<^) = V(qi — q%). In our calculations, we 
consider a one-dimension system and use a Gaussian impurity potential with range a along 
the growth direction 

K(ri-r a ) = V e~^- r ^ 2/a2 . (2.8) 
The Fourier transform of the impurity potential is 

V(qi - q 2 ) = V (27r) 2 a^5(q 1± - d^e^i^nf/^ (2.9) 
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Omitting the perpendicular component k^, q^, we find for the Dyson equation 

G R (k h k[,u;) = G^{k h k\^)+ J d ^Gl(k h q h uo) 

xe -< 9|| - 9 V(gn - <®cr R (j p k> p J). (2.10) 

The above Green's function is for only one impurity at site R and for a single scattering event. 
When all scattering events at that site are included, the Green's function and corresponding 
T-matrix are 

G R (k p k\,u) = G;(k p k' p u) + J d -^Gl{k h q h u) 

xT R (q p q' p uj)G R (q> p k> p u), (2.11) 

xe^^GJ (p||,pj,,w) T' R [p' p ^) . (2.12) 

We assume the impurity density is sufficiently dilute that we can approximate the self energy 
by its impurity-averaged form. The self energy is then 

S r (q\\ , 9|| , = c J dRT R (q\\ , g( , w) , (2. 13) 

where c is the line concentration of impurities. T R differs from T' R of Eq. 1)2.12)1 : we apply 
the self-consistent Born approximatio n 23 ! 27 and we replace G r Ano scattering) in Eq. (J2.11)) by 
G r (including scattering) to obtain T R . We find 

Ta( 9 ,4,J) = e-(»-0V („ -4) + f%f ( 9|| - „) 

--,)"(? ( p , , p j , w ) Ts ( p | h ,J (2. 14) 

The final Dyson equation for the retarded Green's function G r is thus 

CT (k p k' p u)=G r (f> (k p k' p u J ) + J^J^G^(k hn ,Lj) 

xE r (q h g f p u)G r (g' p ^ p u). (2.15) 

Similar equations for the resonant level model (RLM) were obtained by Jauhc 23 . Here 
our approach is rather different. Instead of exploring the relationship of these equations to 
the Boltzmann equation^ 3 , we plan to apply the equation numerically to investigate in detail 
electron-phonon resonances. 
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Eqs. ()2.14j) and ()2.15|) form a closed set of integral equations. In our numerical method, 
we use a uniformly distributed grid for k-space and R-space integration, hence the double 
integral becomes a product of matrices. The sizes of these matrices G r and Tr are typically 
600 x 600 (15 Brillouin zones and 40 points per zone). The R-space integration in Eq. 
(|2.13|) involves about 500 points with a 1000 A grid spacing. To solve the equation set, we 
first use as an initial guess for G r and perform the iteration to obtain a self-consistent 
solution. Sometimes the iterative process does not converge, in which case we adjust the 
initial retarded Green's function, matrix size, or (k,R)-space integral mesh. 

C. Polar optical phonon scattering and G < 

We assume phonon scattering is very weak compare to elastic scattering. In this limit 
the Wannier-Stark level broadening is mainly governed by elastic impurity scattering - 
thus the densities of states are determined by the elastic impurity scattering. We obtain 
the retarded Green's function, which determines the DOS, by the procedure of Sec. Ill Bl 
(neglecting phonon scattering). The approach to G K must be different, however, for this 
function describes the carrier dynamics. When only elastic impurity scattering is considered, 
G K = —G > = —l/2[G r — G a ] = iA/2, and the net current vanishes. In order to obtain a 
non-zero current, phonon scattering must be taken into account in G < . 

We evaluate the lesser Green's function directly from the Keldysh equation^, but with a 
perturbative approach. 

G< = (1 + G r S r )G<(l + S a G a ) + G r S<G a , (2.16) 

where the double momentum labels for each of the Green's functions and self-energies are 
implicit and all the products in the equation imply integration over the adjacent momentum 
label. Here Gq is the lesser Green's function without scattering and is equal to — [GJ — 
G£]/2 = iA(k,ou)/2. The spectral function neglecting scattering, A(k,u) (Eq. (|2.5jl ). has 
several poles of the form S(E — E n ). Scattering makes the poles move away from these 
energies E n ; the first term in Eq. ()2.16|) is proportional to (Gq)~ 1 Gq {Gq)' 1 oc (E—E n ) 2 5(E— 
E n ) near the former singular points. As a consequence, the first term in Eq. ()2.16|) will vanish, 
so Eq. ()2.16|) simplifies to 

G < = G r^< G a ( 2 _ 17 ) 
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Here the self energy is E < 
for impurity scattering is 



+ T, pho . Based on analytic continuation^, the self energy 



^fmp — c I dRT B ^G < T B [. 



(2.18) 



For polar optical phonon scattering, the corresponding self energy is expressed as 
^ho = I [n B G< (w - u q ) + (n B + 1) G< (w + u> q )] , 



(2.19) 



where M q is the electron-phonon coupling matrix element, here treated as a constant, and 
riB = 1/ {exp[hu) q /(kBT)] — 1} is the Bose-Einstein occupation factor for LO phonons. The 
frequency of the LO phonons uj q is also treated as a constant. The temperature we use is 
77K. Earlier wor k 12 i 14 shows higher temperatures will strongly depress the electron-phonon 
resonances. 

The following symmetry holds when there is only elastic scattering: G(k, k',u + ujq) = 



i(k k')—gr 



-G(k,k',uj) or alternately G(k, k', t, t') = G(k,t-t')5[k-k' -^-(t-t')}. We assume 
this symmetry is still present when phonon scattering is included as a perturbation. As a first 
order approximation, we replace G < in Eqs. ()2.18j) and ()2.19j) with \A imp = — | [Gl mp — Gf mp \ . 
The lesser Green's function of Eq. (12.17)1 then becomes 
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The first term has been simplified because iA imp is a solution for the equation G 
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c J dRTftG^npTft G\ . The variables on the right side of Eq. ()2.20|1 are all known, so 
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III. RESULTS 



A. Impurity scattering effects on electron-phonon resonance 

We expect the properties of the electron-phonon resonances to depend significantly on 
the relative magnitudes of the growth axis bandwidth of the lowest conduction miniband 
and the optical phonon energy. If the bandwidth A is very large compared to the optical 
phonon energy, then the electron-phonon resonances can be well described within a Wannier- 
Stark hopping picture, whereas in the other limit the electron-phonon resonances are better 
described within a sequential tunnelling picture. We present here detailed results for the 
drift velocity as a function of electric field for two superlattices, one in each regime. 

The first superlattice consists of 12 monolayers of GaAs and 6 monolayers of AlAs (unit 
cell length 50. 9A). The resulting width of the lowest miniband is 20.3meV. We use the value 
for the GaAs optical phonon energy (Hu = 36meV) for this GaAs/AlAs SL, so the ratio 
of the bandwidth to the optical phonon energy A/hu = 0.56. The other structure is 7 
monolayers of InAs and 12 monolayers of GaSb (unit cell length 57.8A), which has a very 
broad lowest conduction miniband width of 150meV. The optical phonon energy of InAs 
(Two = 30meV) is used for the InAs/GaSb superlattice, and hence A/hu — 5. In both cases, 
the zone center gap between the lowest and second-lowest minibands is larger than 350meV. 
Thus, the occupation of higher minibands is very small and we can neglect conduction from 
higher minibands and interminiband Zener tunneling. For convenience in our comparison 
we approximate the period for the two types of superlattices with the same value, L = 57 A. 

We have numerically calculated the drifty velocity for fields near the electron-phonon 
resonance value neFL = hu, which typically ranges from 0.1 x 10 7 V/m to 1.0 x 10 7 
V/m. Because we treat phonon scattering as a perturbation, the calculated drift velocity is 
directly proportional to the electron-phonon coupling coefficient . We choose the coupling 
constant to make the drift velocity realistic (~ 10 3 m/s). To be consistent, the calculated 
effective relaxation time r e // from phonon scattering should be longer than the relaxation 
time from impurity scattering: r e ff > lOps. 

In Fig. ^a), we show how the density of impurities affects the electron-phonon resonance 
peaks in the narrow bandwith superlattice, and in Fig. E^b) results for the broad bandwidth 
one. The line concentration of impurities ranges from 1/(464A) (defined as c in Fig. ^) 
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to 1/(58A). These values correspond to 1.0 x 10 16 cm~ 3 and 5.1 x 10 18 cm -3 , respectively, 
in 3-dimensions. We see that large impurity densities will produce weaker electron-phonon 
resonance peaks. When the impurity density is 4c, there is only one peak in the velocity- 
field curve for the broad bandwidth SL (Fig. E]^b)). Indeed, the position of the peak is not 
at a resonance, so it's not an electron-phonon resonance at all. Instead, it's the maximum 
of the entire velocity-field curve, which is the starting point of NDC In Fig. ^a), all the 
electron-phonon resonance peaks are located in the NDC region. In Fig. E^b), when the 
impurity scattering is very strong the electric field corresponding to the maximum current 
enters the region of our interest. This peak flattens the electron-phonon resonance peaks. 
The calculated DOS, instead of exhibiting Wannier-Stark levels, is closer to a zero-field DOS. 

In Fig. El the solid lines show the DOS for different impurity concentrations under the 
same electric field F = 0.2 x 10 7 V/m in broad-band superlattice. The DOS for a low impurity 
density (c=l/(464A), thin line in the Fig. |2J) clearly shows the Wannier-Stark levels. But 
the DOS for a higher impurity concentration (c=l/(116A), thick line) tends to be very 
similar to the zero- field result (which is shown as the dashed line). Properly approaching 
the zero-field DOS and the appearance of the maximum of current in the velocity-field curves 
indicate our method is valid for low-field situations. This contrasts with the Wannier-Stark 
hopping method that is often used in high-field transport, which will diverge under low- field 
condition a 12 ! 15 . 

In Fig. |3J we compare the effects of the impurity strength Vq and impurity density c on 
the electron-phonon resonances. We find that the impurity strength can diminish the drift 
velocity considerably, but does not affect the height of the electron-phonon resonances as 
much as the concentration of impurities does. In Fig. 01 we change the impurity size but 
keep Voa fixed to observe changes in the electron-phonon resonances. When the impurity 
size increases, the resonance becomes weaker. 

Figs. H3 and 0] are calculated results for relatively perfect materials. We also model more 
imperfect materials by employing an impurity density of 1/(19A) (corresponding to 1.46 x 
10 20 cm~ 3 in three dimensions), and use a rather small value for the impurity strength Vo- The 
impurity size produces different effects on the electron-phonon resonance peaks for different 
bandwidths. In Fig. |3^a), for the narrow bandwidth SL, impurities of larger size (with Voa 
constant) flatten the electron-phonon resonance peaks. This contrasts with Fig. Efb), for 
the broad bandwidth SL, where larger impurities enhance the electron-phonon resonances 
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stronger. This phenomena can be explained by comparing with the one-dimensional square 
potential barrier scattering problem. The reflection coefficient of a scattering electron is not 
a monotonic function of the barrier size (still keeping V^a constant)^. Instead, there is a 
maximum value for the reflection probability. We take A/2 as the average energy of incoming 
electrons and calculate the impurity size that maximizes the reflection coefficient. For the 
narrow bandwidth SL, a max = 11. 3A, whereas for the broad bandwidth SL, a max = 1.5A. 
Hence impurity scattering generates different trends for narrow and broad bandwidth SLs 
in the region 2A < a < 8A. 

There is another feature distinguishing relatively clean and dirty superlattices. In the 
clean case the electron-phonon resonance in the narrow bandwidth SL is stronger than in 
broad bandwidth one. But in the dirty case, the electron-phonon resonance in the narrow 
bandwidth SL is weaker than in the broad bandwidth one, e.g. when the impurity size is 
a = 8 A. The reason for this is as follows: when V$a is large, all the collisions between 
electrons and impurities are essentially total reflection. In this case, impurity scattering 
in a broad bandwidth SL is stronger than in a narrow bandwidth one because the spatial 
extent of the Wannier-Stark wave function A/eF is proportional to the bandwidth A. Thus 
electrons in a broad bandwidth SL experience more collisions in a Wannier-Stark oscilla- 
tion period. But when V^a is small, the collisions have both substantial transmission and 
reflection. The incoming electrons in a broad bandwidth SL have more energy and a larger 
transmission coefficient than in a narrow bandwidth SL. This makes impurity scattering in 
broad bandwidth SLs weaker than in narrow bandwidth ones, which is the opposite of the 
clean case. 

B. Comparison with earlier results 

We have also reproduced the earlier work of Bryksin and Kleinert^ and find a relation 
between Bryksin's parameter for impurity scattering and ours: hy/U ~ V^y/ac. However, 
this relation has a limited range of applicability because the DOS of a Wannier-Stark level 
generally does not have the half-elliptical shape employed by Bryksin, 

oc l/{2-KhU)Re\/AU - u 2 . (3.1) 



11 



Fig. 13 shows the DOS for different parameter sets and compares them to Bryksin's and 
Wacker' a 14 ! 15 results. When ac is very large, the DOS tends to have a half-elliptical shape. 
From Fig. El we also can see a shift of the DOS. If there is no scattering, the main band 
is centered at E = A/2. Impurity scattering will shift all bands to higher energies. When 
ac is big enough the shape of the DOS is Bryksin-like and the energy shift has the form 
E ce nter = -Reenter + 5E = A/2 + y/nVoac. Bryksin's DOS centers on the zero-energy point 
because they used A — > 0. Wacker's DOS also remains centered on the point E = because 
they chose localized Wanner-Stark states as the basis of the wave function, as have many 
others^. Wacker's DOS also uses the approximation of a constant scattering self-energy, 

rr(F . Jn (A/2eFL) 

( ' k±) -?£-neFL-£ kx + (,r/2)' (3 ' 2) 

where —iT/2 is the self-energy (in Fig. El we adopt T = 2hy/U). The shape of Wannier- 
Stark state using a constant self-energy is not realistic. However, Eq. (|3.2J) does produce the 
appropriate side band. Bryksin's DOS uses the miniband approximation A/2eEL <C 1 and 
neglects the side band, and therefore obtains a field-independent DOS. For the narrow-band 
SL, A =20meV, the central band occupies Jq(A/2cFL) = 90.7% of the entire density when 
a typical electric field F = 0.4 x 10 7 V/m is used. 

In Fig. [?l we reproduce Bryksin's velocity-field curve using a very large value for ac. The 
agreement is very good. When ac is small, the DOS deviates from the half-elliptical shape 
and becomes asymmetric. The asymmetry manifests in a low energy side which has a larger 
density than the high energy side [see Fig. Efa)]. For a realistic superlattice, ac < 1, so the 
half-elliptical shape will typically not exist. 

We have also calculated the approximate momentum relaxation times (r) for different 
impurity parameters and present them in Table 1. The r's are calculated by fitting the 
retarded Green's function \G r (k,t)\ to an exponential decay e _t//T . In Bryksin's limit (using 
the half-elliptical DOS), G r (t) has the form J\(2\/Ut)/VUt, which oscillates in time and 
exhibits a power-law decay. In our calculations, the decay of G r is neither Bryksin-like nor 
exponential. The decay of G r (t) is Bryksin-like when t is short, and the decay is exponential 
when t is long. In the intermediate region, the decay displays a complex shape that depends 
on the specific impurity scattering parameters. 
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IV. DISCUSSION 

The advantages of our computational method are: (i) we treat realistic scattering 
strengths and obtain a realistic DOS; (ii) many common approximations are relaxed, such 
as the small bandwidth (A — > 0), the high-field {^T e ff ^> 1) and the relaxation time ap- 
proximations; (iii) compared to the Wannier-Stark hopping method that was mainly used 
under high- field conditions, our method is applicable to both high and low field transport; 
and (iv) we do not use the generalized Kadanoff-Baym ansatz 18 to solve the Dyson equation. 

In our calculations, we used the approximation that phonon scattering is weak compared 
to impurity scattering. It is straightforward in principle to eliminate this approximation 
in the non-equilibrium Green's function method by adding self-energy terms for phonon 
scattering TT h ^- to Ea. fHoT) . 



If we ignore the last term in Eq. (|4.1j) . which corresponds to a low carrier concentration 
approximation, then Eq. (J2.15|) is a closed equation for the retarded Green's function. Based 
on the retarded Green's function, we can use iteration or some other method to solve the 
Quantum Kinetic Equation (j2.17|) to obtain the lesser Green's function. If the low carrier 
concentration approximation is relaxed, Eq. (j2.15|) and Eq. (j2.17|) are coupled and more 
iterations are needed. In the both cases, the computational effort is much greater than the 
present approach. 
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Table 1. Calculated approximate relaxation times for different impurity scattering 
strengths. 



Impurity parameters 


Approximate relaxation time(ps) 


a(A) 


VO(eV) 


c(l/A) 


A = 20meV 


A = 150meV 


20 


0.1 


1/464 


8.95 


1.73 


20 


0.1 


2/464 


3.09 


0.377 


20 


0.1 


4/464 


0.905 


0.116 


20 


0.1 


8/464 


0.470 




100 


0.02 


1/464 


3.32 


1.18 


100 


0.04 


1/464 


2.12 


0.683 


100 


0.02 


2/464 


0.940 


0.321 
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Figure captions 



Figure 1: Drift velocity versus field for different impurity concentration (a) in a narrow 
band superlattice (A = 20meV) and (b) in a broad band superlattice (A = 150meV). The 
unit of impurity concentration c = 1/(464A). Inset of (a) shows a clear view of the velocity- 
field curves for the high impurity densities 4c and 8c because the drift velocity is very small. 
The impurity scattering strength Vq = O.leV, and the impurity size a = 

Figure 2: The calculated DOS for different impurity concentrations for the broad band- 
width superlattice (A = 150m V). c = 1/(464A) and the electric field F is in units of 10 7 V/m. 
The zero-field DOS is shown for comparison (dashed line). 

Figure 3: Different effects of impurity strength and concentration on electron-phonon 
resonance peaks in a (a) narrow bandwidth and (b) broad bandwidth superlattice. The 
parameters used are c = 1/(464A), V = 0.02eV, and a = lOOA. 

Figure 4: velocity-field curves for different impurity sizes with a fixed VqCl (= 2eV-A) in 
a (a) narrow bandwidth and (b) broad bandwidth superlattice. The impurity concentration 
c= 1/(464 A). 

Figure 5: The drift velocity-field relations for different impurity sizes while Vqcl keeps 
constant (= 0.12eV- A) in (a) narrow bandwidth and (b) broad bandwidth superlattices. 
The impurity density c = 1/(19A). 

Figure 6: (a) calculated DOS for different ac, for Voy/ac = const. (1.7meV), (b) Bryksin- 
type DOS (hVU = 1.7meV); and Wacker-type DOS, with a constant self-energy T/2 = 
hVU = 1.7meV. 

Figure 7: Reproduction of Bryksin's results for two bandwidths (A = 20meV, 36meV). 
Bryksin's parameter h\/U = 1.7meV. Our parameters are V = 0.425meV and ac = 16. 
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